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Abstract 

Stability analysis of agegraphic dark energy in Brans-Dicke theory is presented in this paper. We 
constrain the model parameters with the observational data and thus the results become broadly 
consistent with those expected from experiment. Stability analysis of the model without best fitting 
shows that universe may begin from an unstable state passing a saddle point and finally become 
stable in future. However, with the best fitted model, There is no saddle intermediate state. The 
agegraphic dark energy in the model by itself exhibits a phantom behavior. However, contribution 
of cold dark matter on the effective energy density modifies the state of teh universe from phantom 
phase to quintessence one. The statefinder diagnosis also indicates that the universe leaves an 
unstable state in the past, passes the LCDM state and finally approaches the sable state in future. 
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1. INTRODUCTION 



Cosmological probes such as Cosmic Microwave Background (CMB) B], Supernova 
type la (SNIa)Q, Weak LensingQ], Baryon Acoustic Oscillations (BAO) 6|, 2dF Galaxy 
Redshift Survey (2dFGRS) j^j at low redshift and DEEP2 redshift survey 9] at high red- 
shift, have given us cross-checked data to determine cosmological parameters with high 
precision. These parameters indicate that our universe is nearly spatially flat, homogeneous 



ker (FRW) with zero curva- 
In addition, from AC DM 



and isotropic at large scale, i.e. a Friedmann-Robertson-Wai 
ture, and has entered an accelerating epoch at about z ~ 0.46 
model, the universe consists of 0.046 baryonic matter, 0.228 non-relativistic unknown matter, 
called as cold dark matter (CDM), and a significant amount of 0.726 smoothly distributed 
dominant dark energy (DE) Q]. 

The equation of state (EoS) of DE, is the main parameter which determines the grav- 
itational effect on the dynamics of the universe, and can be measured from observations 
without need to have a definite model of DE. Strong evidences suggest that the EoS of DE 
lies in a narrow range around w ~ — 1 with a smooth behavior jg], 10 1. We can classify the 
EoS parameter of DE with respect to the barrier w = —1, namely the phantom divide line 
(PDL ) 11] . That is, DE with the EoS parameter of w = — 1 is named for the cosmological 
constant, A, with a constant energy density. The case with dynamical EoS parameter o 



w > — 1 and w < — 1 are respectively referred to quintessence and phantom cosmology 12 1 

A fundamentally important problem in modern cosmology is the current state of the 
universe to be in phantom or quintessence phase? A transition from phantom to quintessence 
phase and reverse plus variation of EoS parameter being close to -1 imply that universe 
may evolve into a finite-time future singularity in future 15 j isl - Different cosmolo gica l 



models predict different scenarios based on their own assumptions and criterions [18] 19] 20]. 
Among them, the inclusion of CDM in the formalism, and its coupling to DE and or geometry 
to explain phantom transition, varying EoS parameter and future singularity are favored by 
many authors 49J 22]. 

Although, the issue of DE, its energy density and EoS parameter is still an unsolved 
problem in classical gravity, it may be in the context of quantum gravity that we achieve 
a more inclusive insight to its properties 27]. The holographic dark energy (HDE) model 



2 



is an alternative to formulate DE within the framework of quantum gravity |28[, 29|. The 
holographic principle expresses that the number of degrees of freedom describing the physics 
inside a volume (including gravity) is confined by the area of the boundary which encloses 



the volume and therefore related to the entropy scales with the boundary [30|. The entropy 
scales like the area rather than the volume, and thus the fundamental degrees of freedom 
rendering the system are characterized by an effective quantum field theory in a box of size 



30]. 



L with one fewer space dimensions and with planck-scale UV cut-off A 

On the other hand, a viable alternative approach to dark energy that naturally gives rise 
to late time accelerating solutions as well as early time inflation is extended modified gravity 
(For a review see ). Among extended theories of gravity, scalar tensor theories are the best 
motivated and most promising candidates to general relativity They are compatible with 
the observational findings without the need of additional components and may resolve the 
coincidence problem. The birth of scalar tensor theories is marked by the work of Jordan in 
which a scalar field coupled to both curvature and matter Lagrangian. Later, Brans-Dick 
demanded the decoupling of matter Lagrangian from scalar field to guarantee the validity 
of the weak equivalence principle. Though, the "decoupling of the scalar" (called as Brans- 
Dicke problem) is the most serious shortcoming of this approach, it has been argued that 
the cosmic acceleration can be generated by either assuming a time dependent Brans-Dicke 



23||2^|25||46|. 



parameter out, or adding a potential term to the lagrangian 

Alternative, with regards to HDE, the scalar-tensor theories have been widely used to 
explain the late time acceleration of the universe js3|-62]. In particular, in Brans-Dicke 
(BD) theory 69] 84J 85] 23] 83], where the gravitational constant varying as inverse of a 
time dependent scalar field, cosmic acceleration is predicted. 

While the successful HDE model explains the observational data, more recently, a new 
dark energy model, dubbed agegraphic dark energy (ADE) model, has been proposed by 



Cai 



38 



]. The ADE is also related to the holographic principle of quantum gravity and 
considers the uncertainty relation of quantum mechanics together with the gravitational 
effect in general relativity. 

In this work we present the ADE in Brans-Dicke theory after performing stability anal- 
ysis and constraining the model parameters with observational data. As noted by 2Q|, the 
acceptable cosmic acceleration in BD theory can be realised only if the BD parameter is 
time dependent, or a potential term included in the action. However, in this work, we 
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reproduce current universe acceleration in the context of BD theory without adding po- 
tential function or assuming time dependent BD parameter. Instead, with the agegraphic 
description of dark energy, the best fitted effective EoS parameter exhibits current cosmic 
acceleration in quintessence phase and matter dominated phase in the far past. Together 



with statefinder parameters [7l| 72J- 75| 76|, we also discuss the dynamics of the model in 
different cosmological epochs. 

2. THE MODEL 

The Brans-Dicke action is given by, 

S = J d 4 x^ (f-R - i<rv^v„0 + L A ?j , (1) 

where u is a dimensionless coupling constant which determines the coupling between gravity 
and BD scaler field, R is Ricci scaler, is the BD scalar field and Lm Lagrangian of 

matter field. We assume the metric is in the form of, 

ds * = + tffa y) f _^!_ + t*<kA , (2) 

where space time is homogeneous and isotropic (FRW universe). In the metric, k is curvature 
parameter with k = —1,0, 1 corresponding to open, flat, and closed universes, respectively. 
The BD scaler field <fi and the scale factor a, are functions of t only. Variation of action (jTJ) 
with respect to metric ([2]) yields the following field equation, 

^ (2* + J/ 2 + A) - -HH> -±M>- 

^fl + -V=PA, (4) 



2 V u 

4> + 3H<) ) -^(- + HZ + ^) ( f ) = 0, (5) 
2uj\a or / 

where H = -. Moreover, the conservation equations for dark energy and matter field in the 
universe are respectively, 

p A + 3Hp A (l + w A ) = 0, (6) 
p m + 3Hp m = 0, (7) 



Next, we apply the ADE model in Brans-Dicke theory. The ADE model with dark energy 
density is given by 



Pa 



3n 2 M 2 



J^2 

where T is age of the universe and given by 



In the framework of Brans-Dicke cosmology, we write the ADE of quantum fluctuations in 
the universe as 

P» = (10) 

In the next section we constrain the model with the data from observational probe, and 
study the stability of the dynamical system. 



3. STABILITY ANALYSIS AND OBSERVATIONAL CONSTRAINTS 



The structure of the dynamical system can be studied via phase plane analysis. By 
introducing the following dimensionless variables, 



a 



4up n 



k 



a- 



n 



(11) 



30 2 # 2 ' # 2 a 2 ' 11 H 2 T 2 

We assume that the BD scalar field is in power law of the scale factor a(t), i.e. (j)(t) oc a a . 
Therefore, the field equations in terms of new dynamical variables become, 



-fi m (3 + 2« + 2— ) 



-2Q 



D 



n 



H 2 ' 



where 



H 

H 2 



2(3a + 3) 

+n k - i + 30 A (— i 



(2a + 1) 
2a 



■ 2a(au 

2y/Q^- 



3 
2 



3 3n 



3 
2 



(12) 
(13) 
(14) 



(15) 
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By using the Friedman constraint in terms of the new dynamical variables: 

au s 



Q fc = Vt m + VL A - 2a(l 



1 



the equations (jl2l!14p reduce to, 

SI' = Sir 



2a 



(3a + 3) 



((2a + l) 2 + a 



3, , «u, 

+ 2a(au - -) + n A - 2a(l - — ) 



- 2 + 3Sl A (-l 



2a 2^/ST, 



3n 



(16) 



(17) 



si, 



-2x/n f 



■n 



(3a + 3) 



(2a + r 



+ 2a(au --)- 2a(l - — ) - 2 + Q m + Q A 



3Q A (-1 



2a 
~3~ 



2v^ 

3n 



(18) 



where prime means derivative with respect to ln(a). In stability formalism, by simultane- 
ously solving Q' m = and Q' A = 0, we find the fixed points (critical points). However, in 
this work we first constrain the model with observational data for distance modulus using 
X 2 method. The model parameters we best fit are a, and the ADE, BD and Hubble pa- 
rameters, n, u and ho respectively. The observational data we use is the most recent SNe 
la data given by 557 data points. The chi 2 method will be employed to best fit the model 
parameters with the observational data. Figs 1 and 2 show the best fitted 2-dim likelihood 
and confidence level for the parameters a, ho, n and oj. 
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Fig. l:The best-fitted 2-dim likelihood and confidence level for a and h 




Fig. 2:The best-fitted 2-dim likelihood and confidence level for n and u 

We the values of the best fitted parameters we study the stability of the model. We 
solve the dynamical system and obtain the critical points. From eigenvalues, the stability 
of the critical points are given in Table I. From Table I, there are five critical points where 



TABLE I: Best fitted critical points 



point 




eigenvalue 


stability 


Pi 


(4.7,0) 


(2.8,4.4) 


unstable 


P2 


(0,0.6) 


(-1.1,-3.9) 


stable 


P3 


(0,0) 


(-2.8,1.6) 


saddle 


Pi 


(0,27.8) 


(1.9,-0.9) 


saddle 


P 5 


(-38.8,43.6) 


(3 x 10^7,2.8 


) unstable 



only one of them is stable. In the 2-dim phase space ( Fig 3), the best fitted trajectory 
leaves unstable critical point Pi, towards the stable critical point P% is shown. Note that 
the critical points are presenting the state of the universe. By small perturbation of the 
dynamical variables, the green color best fitted trajectory in the phase plane illustrates that 
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the dynamical universe begins from an unstable state Pi in the past in the phase plane and 
approaches the stable state P2 in future. 
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Fig. 3:The 2-dim phase plane corresponding to 
the critical point 

At this point, from the graph of phase plane not enough physical information can be 
found. Our original model is presented in term of new dimensionless variables with not 
enough physical justification. In addition, the phase space shows the behavior of the system 
in terms of these new variables. One needs to represent the physical quantities in terms 
of these variables or alternatively switch back into the original variables to understand the 
cosmological consequences. More specifically, in the following we discuss the cosmologi- 
cal parameters such as EoS parameter and statefinders that are among the most practical 
parameters to discriminate cosmological models. The effective EoS parameter, w e ff, age- 
graphic EoS parameter, are in particular given in terms of new dimensionless variables 
as 



1 



W eff 



3a + 3 
-2a(l - 



((2a + l) 2 + 2a(au - - 



-2 + 3fi A (-l 



2' 
2a 



f2 m -|- 



+ 



2JTt, 



3n 



(19) 



W A 



2a 



2Vfi A 
3n 



(20) 



Also, the statefinder parameters in terms of deceleration parameter q and are defined as 
r = § - 3g - 2 and s = 

In Table II, the values of the effective EoS and deceleration parameters in addition to 
statefinder parameters {r, s} are shown at the critical points and also at current epoch. 



S 



TABLE II: The values of cosmological parameters 



parameter w e ff q r s 



current value -0.8 -0.695 0.3 0.17 
unstable Point 0.467 1.2 4.08 -1.34 
stable Point -0.82 -0.732 0.35 0.14 



The effective EoS parameters is shown in Fig. 4)top. The parameter begins to change 
from unstable critical point in the past moving towards the stable critical point in the future. 
From the graph, its current value is —0.8. A comparison of the effective EoS parameter and 
ADE EoS parameter is given in Fig. 4)below. From the graph we see that w\ is always less 
than -1 and thus represents the phantom regime. The contribution of CDM in the universe 
causes the effective EoS parameter levels up from phantom phase to quintessence phase. 
Note that these trajectories are for the best fitted parameters with the observational data 
and therefore are physically significant. 
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Fig. 4: The graph of effective EoS parameter w e ff, as functions of redshift 

Fig 5 shows the trajectories of the statefinder diagrams {r, q}, {s,q} and {r, s}. From 
the graph it can be seen that the best-fitted trajectory passes LCDM state with r,s = 1,0 
sometimes in the past. The current value of the best fitted trajectory and its location with 
respect to the LCDM state can also be observed in the r, s diagram. 
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Fig. 5: The plot of statefinder parameters {r, q}, {s, q} and {r, s} 

Moreover, Fig 6 shows the corresponding dynamical behavior of the satefinder r and s as 
a function of N = ln(a). 



r 




-10 -5 5 10 -10 -5 5 10 



N N 

Fig. 6: The plot of statefmder parameters r and s as function of N = ln(a). 



4. SUMMARY AND DISCUSSION 



The paper is designed to investigate the agegraphic dark energy model in Brans-Dicke 
theory. After the equations have been cast in the appropriate form, from stability analysis 
and best fitting parameters, one stable and four unstable critical points are obtained. The 
best fitted trajectory shows that the universe moves from unstable state in the past to the 
attractor state in future. Testing the model against EoS parameter for agegraphic dark 
energy implies that universe is always in phantom era. However, if contribution of CDM 
into the total energy to be taken into account, then the effective EoS parameter reveals that 
the universe at higher redshifts is in matter dominated epoch and just recently at about 
z = 0.7, it enters quintessence phase. The statefmder parameters also show that universe 
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has passed LCDM state in the past and approaches stable state in future. 



[1] Dunkley J. et al., (2009), Astrophys. J. 701, 1804 ; Gold B. et al, (2009), Astrophys. J. Suppl. 
180, 265; Hill R.S. et al, (2009), Astrophys. J. Suppl. 180, 246; Hinshaw G. et al, (2009), 
Astrophys. J. Suppl. 180, 225 ; Nolta M.R. et al, (2009), Astrophys. J. Suppl. 180, 296. 

[2] Komatsu E. et al, (2009), Astrophys. J. Suppl. 180, 330. 

[3] Knop R.A. et al., (2003), Astrophys. J. 598 102 ; Garnavish P.M. et al., (1998), Astrophys. 
J. 493,L53; Perlmutter S. et al., (1997), Astrophys. J. 483,565 ; Frieman J.A. et al., (2008), 
Astrophys. J. 135,338 ; Sako M. et al, (2008), Astrophys. J. 135,348. 

[8] Riess A.G. et al, (2004), Astrophys. J. 607,665. 

[5] Leauthand A. et al., (2010), Astrophys. J. 709,97; Kubo J.M. et al., (2009), Astrophys. J. 702, 

L110; Sato M. et al., (2009), Astrophys. J. 701,945; 
[6] Parkinson D. et al., (2010), Mon. Not. Roy. Astron. Soc. 401,2169; Percival W. et al., (2010), 

Mon. Not. Roy. Astron. Soc. 401,2148; Wang X. et al., (2009), Mon. Not. Roy. Astron. Soc. 

394,1775; Benitez N. et al., (2009), Astrophys. J. 691,241. 
[7] Cole S. et al., (2005), Mon. Not. Roy. Astron. Soc. 362, 505; Croton D.J. et al., (2005), Mon. 

Not. Roy. Astron. Soc. 356, 1155; Wild V. et al., (2005), Mon. Not. Roy. Astron. Soc. 356,247. 
[8] Riess A.G. et al., (2004), Astrophys. J. 607,665. 

[9] Yan, R. et al., (2009), Mon. Not. Roy. Astron. Soc. 398,735; Sawicki M. et al, (2008), Astro- 
phys. J. Astrophys. J. 687,884; Schiavon R.P. et al., (2006), Astrophys. J. 651,L93. 

[10] Amanullah R. et al.. larXiv:1004.171H yl. to appear in (2010), Astrophys. J. 

[11] Cai Y.-F. et al., (2010), Phys. Rept., in press, doi:10.1016/j.physrep. 2010.04.001. 

[12] Gonzalez T., Leon G, and Quiros I., (2006), Classical Quantum Gravity 23,165; Johri V.B., 
(2002), Classical Quantum Gravity 19,5959; Brax P. and Martin J., (2000), Phy. Rev. D 
61,103502. 

[13] Carroll S.M., De Felice A., and Trodden M., (2005), Phys. Rev. D 71,023525; Stefancic H., 

(2004) Phys. Lett. B 586,5. 
[14] Caldwell R.R., Kamionkowski M., and Weinberg N.N., (2003), Phys. Rev. Lett. 91, 071301 ; 

Caldwell R.R., (2002), Phys. Lett. B 545,23. 
[15] S. Nojiri, S. D. Odintsov, Phys. Lett. B562:147-152,(2003) 

11 



[16] E. Elizalde, S. Nojiri, S. D. Odintsov, Phys. Rev. D70:043539,(2004) 

[17] S. Capozziello, S. Nojiri, S.D. Odintsov,Phys.Lett.B632:597-604,(2006) 

[18] S. Nojiri, Sergei D. Odintsov, Phys. Rev. D72:023003,(2005) 

[19] S. Nojiri, S. D. Odintsov, Phys. Rev. D70:103522,(2004) 

[20] S. Nojiri, S. D. Odintsov, Phys. Lett. B639:144-150,(2006) 

[49] S. Nojiri, S. D. Odintsov, S. Tsujikawa, Phys. Rev. D71:063004,(2005) 

[22] Ito, S. Nojiri, S. D. Odintsov. larXiv:1111.5389V l 

[23] E. Elizalde, S. Nojiri, S. D. Odintsov, D. S. Gmez, V. Faraoni, Phys. Rev. D 77:106005,(2008) 

[24] S. Nojiri, S. D. Odintsov, Phys. Rept. 505:59-144,(2011) 

[25] S. Nojiri, S.D. Odintsov, ECONF C0602061:06,2006; Int.J.Geom.Meth.Mod.Phys.4:115- 
146,(2007) 

[46] S. Capozziello, S. Nojiri, S.D. Odintsov, A. Troisi, Phys. Lett. B639:135-143,(2006) 

[27] Cohen A. G., Kaplan D. B. and Nelson A. E., (1999), Phys. Rev. Lett. 82,4971 . 

[28] Hsu S. D. H., (2004), Phys. Lett. B 594,13. 
Li M., (2004), Phys. Lett. B 603,1. 



[30] Hooft G. t , |gr-qc/9310026| ; Susskind L., (1995), J. Math. Phys, 36,6377. 

[31] Cai Y.-F. and Wang J., (2008), Class. Quant. Grav. 25,165014; Guo Z.-K. et al., (2005), Phys. 

Lett. B 608,177; Feng B., Wang X.L., and Zhang X.M., (2005), Phys. Lett. B 607,35 . 
[32] Xia J.Q. et al., (2009), Int. J. Mod. Phys. D 17, 1229 ; Kunz M. and Sapone D., (2006), Phys. 

Rev. D74, 123503; Caldwell R.R. and Doran M., (2005), Phys. Rev. D. 71,23515 . 
[33] Copeland E.J., Sami M., and Tsujikawa S., (2006), Int. J. Mod. Phys. D 15, 1753 ; Peebles 

P.J.E. and Ratra B., (2003), Rev. Mod. Phys. 75,559. 
[34] Qiang, L.Y.Ma, and Yu D., (2005), phys. Rev. D 71,061501. 

[35] Huang Q. G. and Li M., (2004), J. Cosmol. Astropart. Phys. 0408, 013 ;(2005), J. Cosmol. As- 
tropart. Phys. 0503, 001; Setare M. R., Zhang J. and Zhang X., (2007), J. Cosmol. Astropart. 
Phys. 0703,007. 

[36] Wang B., Gong Y. G. and Abdalla E., (2005), Phys. Lett. B 624, 141; Wang B., Lin C. Y. 

and Abdalla E., (2006), Phys. Lett. B 637,357. 
[37] Zhang X. and Wu F. Q., (2005), Phys. Rev. D 72, 043524; (2007), Phys. Rev. D 76, 023502; 

Chang Z., Wu F. Q., and Zhang X., (2006), Phys. Lett. B 633,14. 
[38] Cai R. G., (2007), Phys. Lett. B 657,228 . 



12 



[39] Overdium J. M. and Wesson P. s., (1997), phys. Rep. 283,303. 
[83] Sahoo B.K. and Singh L.P., (2002), Mod. Phys. Lett. A 17,2409. 
[41] Boisseau B. et al., (2000), Phys. Rev. Lett. 85,2236. 
[42] Gannouji R. et al., (2006), J. Cosmol. Astropart. Phys. 0609,016. 
[43] Sahoo B.K. and Singh L.P., (2003), Mod. Phys. Lett. A 18, 2725. 
[44] Sahoo B.K. and Singh L.P., (2004), Mod. Phys. Lett. A 19, 1745. 

[45] Sadeghi J., Setare M.R., Banijamali A. and Milani F., (2009), Phys. Rev. D 79, 123003. 
[46] Capozziello S., Carloni S. and Troisi A., (2003), Recent Res. Dev. Astron. Astrophys. 1 625. 
[47] Nojiri S. and Odintsov S.D., (2003), Phys. Rev. D 68 123512; (2003), Phys. Lett. B 576, 5; 

Farajollahi H., Salehi A., Tayebi F., Ravanpak A., (2011), J. Cosmol. Astropart. Phys. 05, 

017. 

[48] Faraoni V., (2007), Phys. Rev. D 75 067302; Briscese F., Elizalde E., Nojiri S. and Odintsov 

S.D., (2007), Phys. Lett. B646, 105. 
[49] Nojiri S. and Odintsov S.D. , (2004), Gen. Rel. Grav. 36, 1765 ; (2004), Phys. Lett. B 599, 137 

; Farajollahi H., Milani F., (2010), Mod. Phys. Lett. A 25:2349-2362; Farajollahi H., Setare 

M. R., Milani F. and Tayebi F., (2011), Gen.Rel.Grav.43:1657-1669. 
[50] Setare M. R., Jamil M., (2010), Phys. Lett. B 690 1-4 ; Ito Y., Nojir S. i, (2009), 

Phys.Rev.D79: 103008. 
[51] Mota D.F., Shaw D.J., (2007), Phys. Rev. D 75, 063501. 
[52] Dimopoulos K., Axenides M., (2005), J. Cosmol. Astropart. Phys. 0506:008. 
[53] Damour T., Gibbons G. W. and Gundlach C, (1990), Phys. Rev. Lett, 64, 123.; Farajollahi 

H., Mohamadi N., (2010), Int. J. Theor. Phys.49:72-78; Farajollahi H., Mohamadi N. , Amiri 

H., (2010), Mod. Phys. Let. A, Vol. 25, No. 30 25792589. 
[54] Carroll S. M., (1998), Phys. Rev. Lett. 81 3067. 

[55] Carroll S. M., Press W. H. and Turner E. L., (1992), Ann. Rev. Astron. Astrophys, 30, 499. 



[56] Biswas T. and Mazumdar A., arXiv:hep-th/0408026 



[57] Biswas T., Brandenberger R., Mazumdar A. and Multamaki T., (2006), Phys. Rev. D74 , 
063501. 

[71] Sahni V., Saini T. D., Starobinsky A. A., Alam U., (2003), JETPLett. 77:201-206. 

[76] Linder E.V., (2008), Rep. Prog. Phys. 71, 056901. 

[60] Khoury J. and Weltman A., (2004), Phys. Rev. Lett. 93, 171104. 

13 



[61] Brax P., de Brack C van , Davis A. C, Khoury J., Weltman A., (2004), Phys.Rev.D70:123518 

; Khoury J., Weltman A., (2004), Phys. Rev. D 69, 044026. 
[62] Farajollahi H., Salehi A., (2010), Int. J. Mod. Phy. D, 19(5), 113. 
[63] Bouhmadi-Lopez M., (2008), Nucl. Phys. B. 797, 78-92 . 

[64] Agarwal N., Bean R., Khoury J. and Trodden M., (2010), Phys. Rev. D. 81, 084020 . 

[65] Lee H. M. and Tasinato G., (2004), J. Cosmol. Astropart. Phys. 0404, 009 . 

[66] Karasik D. and Davidson A., (2004), Class. Quant. Grav. 21, 1295-1302 . 

[67] Aoyanagi K. and Maeda K., (2006), J. Cosmol. Astropart. Phys. 0603, 012 . 

[68] Shtanov Y., Viznyuk A. and Sahni V., (2007), Class. Quant. Grav. 24, 6159-6190 . 

[69] Brans C. H. and Dicke R. H., (1961), Phys. Rev. 124, 925. 

[70] N. Banerjee and D. Pavon, Phys. Rev. D 63, 043504 (2001). 

[71] Sahni V., Saini T. D., Starobinsky A. A., Alam U., (2003), JETPLett. 77:201-206. 
[72] Alam U., Sahni V., Saini T. D. and Starobinsky A. A., (2003), Mon. Not. Roy. Astron. Soc. 
344, 1057. 

[73] Zimdahl W., Pavon D., (2004), Gen. Rel. Grav. 36, 1483. 
[74] Yi Z. L. and Zhang T. J., (2007), Phys. Rev. D75, 083515. 

[75] Farajollahi H., Salehi A., (2011), Phys.Rev.D83: 124042, ; Farajollahi H., Salehi A. , Tayebi 
F., Ravanpak A., (2011), JCAP 05 017 ; Farajollahi Hossein, Salehi Amin, (2010), JCAP 
1011:006,; Farajollahi H., Salehi A., Tayebi F., will be published in Astrophys. Space Sci 

[76] Linder E.V., (2008), Rep. Prog. Phys. 71 056901. 

[77] La D., Steinhardt P. J. and Bertschinger E. W., (1989), Phys. Lett. B 231, 231. 

[78] Farajollahi H. and Mohamadi, N., (2010), Int. J. Theor. Phy, 49, 1, 72-78. 

[79] Bento M. C, Bertolami O., and Sa P. M., (1991) Phys. Lett. B 262, 11. 

[80] Bento S. J. and Eardly D. M., (1995), Ann. Phys. (NY.) 241, 128. 

[81] Barrow J. D. et al., (1993), Phys. Rev. D 48, 3630. 

[82] Barrow J. D. et al., (1992), Mod. Phys. Lett. A 7, 911. 

[83] Sahoo B. K. and Singh L. P., (2003), Mod. Phys. Lett. A 18, 2725. 

[84] Mathiazhagan C. and Johri V. B., (1984), Class. Quantum Grav. 1, L29. 

[85] La D. and Steinhardt P. J., (1989), Phys. Rev. Lett 62, 376. 

[86] Perlmutter S. et al., (1999), Astrophys.J. 517, 565. 

[87] Garnavich P. M. et al., (1998), Astrophys. J. 509, 74. 

11 



[88] Bertotti B., less L. and Tortora P., (2004), Nature 425, 374. 



15 



